Interaction of Short-Chain PFAS with Polycationic Gels: How Much Fluorination is Necessary for Efficient Adsorption?

The short-chain per- and polyfluorinated alkyl substances (PFAS), introduced to replace the legacy PFAS compounds, turned out to be as toxic and harmful as their longer-chain predecessors and even harder to sequester from contaminated water sources. In this work, molecular dynamics (MD) simulations are employed to investigate the adsorption mechanism of GenX, a representative compound for short-chain PFAS, on a polycationic hydrogel with various extents of fluorination in its backbone and cross-linkers. Simulations indicate that the presence of fluorinated segments next to cationic groups in the polymer gel structure provides the most efficient environment for GenX adsorption. The combination of electrostatic and hydrophobic interactions offered by the cationic-fluorophilic segments amplifies the binding of GenX molecules compared to polymer segments with nonfluorinated cationic or noncationic fluorinated segments. Moreover, such a gel demonstrates high selectivity toward GenX against its hydrogenated analogue.

Green Polymer Electrolytes Based on Polycaprolactones for Solid-State High-Voltage Lithium Metal Batteries.

Solid polymer electrolytes (SPEs) have attracted considerable attention for high energy solid-state lithium metal batteries (LMBs). In this work, potentially ecofriendly, solid-state poly(ε-caprolactone) (PCL)-based star polymer electrolytes with cross-linked structures (xBt-PCL) are introduced that robustly cycle against LiNi0.6 Mn0.2 Co0.2 O2 (NMC622) composite cathodes, affording long-term stability even at higher current densities. Their superior features allow for sufficient suppression of dendritic lithium deposits, as monitored by 7 Li solid-state NMR. Advantageous electrolyte|electrode interfacial properties derived from cathode impregnation with 1.5 wt% PCL enable decent cell performance until up to 500 cycles at rates of 1C (60 °C), illustrating the high potential of PCL-based SPEs for application in high-voltage LMBs.

GenX in water: Interactions and self-assembly.

2,3,3,3-tetrafluoro-2-(heptafluoropropoxy) propanoate, a.k.a. "GenX", is a surfactant introduced as a safer alternative to replace perfluorooctanoate (PFOA) in the manufacturing of fluorinated polymers, however, GenX is shown to cause adverse health effects similar to, or even worse than, those of the legacy PFOA. With an overarching goal to understand the behavior of GenX molecules in aqueous media, we report here on GenX micelle formation and structure in aqueous solutions, on the basis of results obtained from a combination of experimental techniques such as surface tension, fluorescence, viscosity, and small-angle neutron scattering (SANS), and molecular dynamics (MD) simulations. To our best knowledge, this is the first report on GenX micelles. The critical micelle concentration (CMC) of GenX ammonium salt in water is 175 mM. GenX forms small micelles with association number 6-8 and 10 Å radius. GenX molecules prefer to align along the micelle surface, and the ether oxygen of GenX has very little interaction with and exposure to water. Information on the surfactant and interfacial properties of GenX is crucial, since such properties are manifestations of interactions between GenX molecules and between GenX and water molecules and, in turn, the amphiphilic character of GenX dictates its fate and transport in the aqueous environment, its interactions with various biomolecules, and its binding to adsorbent materials.

Structure and composition of mixed micelles formed by nonionic block copolymers and ionic surfactants in water determined by small-angle neutron scattering with contrast variation.

HYPOTHESIS: Complex fluids comprising polymers and surfactants exhibit interesting properties which depend on the overall composition and solvent quality. The ultimate determinants of the macroscopic properties are the nano-scale association domains. Hence it is important to ascertain the structure and composition of the domains, and how they respond to the overall composition. EXPERIMENTS: The structure and composition of mixed micelles formed in aqueous solution between poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO) block copolymers (Pluronics or Poloxamers) and the ionic surfactant sodium dodecylsulfate (SDS) are determined from an analysis of small-angle neutron scattering (SANS) intensity data obtained at different contrasts. Different polymers and concentrations have been probed. FINDINGS: The SDS + Pluronic mixed micelles include polymer and some water in the micelle core that is formed primarily by alkyl chains. This is different than what was previously reported, but is consistent with a variety of experimental observations. This is the first report on the structure of SDS + Pluronic P123 (EO19PO69EO19) assemblies. The effects on the mixed micelle structure and composition of the surfactant concentration and the polymer hydrophobicity are discussed here in the context of interactions between the different components.

Controlling the self-assembly of perfluorinated surfactants in aqueous environments.

Surface active per- and polyfluoroalkyl substances (PFAS) released in the environment generate great concern in the US and worldwide. The sequestration of PFAS amphiphiles from aqueous media can be limited by their strong tendency to form micelles that plug the pores in the adsorbent material, rendering most of the active surface inaccessible. A joint experimental and simulation approach has been used to investigate the structure of perfluorooctanoate ammonium (PFOA) micelles in aqueous solutions, focusing on the understanding of ethanol addition on PFOA micelle formation and structure. Structurally compact and slightly ellipsoidal in shape, PFOA micelles in pure water become more diffuse with increasing ethanol content, and break into smaller PFOA clusters in aqueous solutions with high ethanol concentration. A transition from a co-surfactant to a co-solvent behavior with the increase of ethanol concentration has been observed by both experiments and simulations, while the latter also provide insight on how to achieve co-solvent conditions with other additives. An improved understanding of how to modulate PFAS surfactant self-assembly in water can inform the fate and transport of PFAS in the environment and the PFAS sequestration from aqueous media.

Structure and Interactions in Perfluorooctanoate Micellar Solutions Revealed by Small-Angle Neutron Scattering and Molecular Dynamics Simulations Studies: Effect of Urea.

The self-assembly of surfactants in aqueous solution can be modulated by the presence of additives including urea, which is a well-known protein denaturant and also present in physiological fluids and agricultural runoff. This study addresses the effects of urea on the structure of micelles formed in water by the fluorinated surfactant perfluoro-n-octanoic acid ammonium salt (PFOA). Analysis of small-angle neutron scattering (SANS) experiments and atomistic molecular dynamics (MD) simulations provide consensus strong evidence for the direct mechanism of urea action on micellization: urea helps solvate the hydrophobic micelle core by localizing at the surface of the core in the place of some water molecules. Consequently, urea decreases electrostatic interactions at the micelle shell, changes the micelle shape from prolate ellipsoid to sphere, and decreases the number of surfactant molecules associating in a micelle. These findings inform the interactions and behavior of surface active per- and polyfluoroalkyl substances (PFAS) released in the aqueous environment and biota.

Small Groups, Big Impact: Eliminating Li+ Traps in Single-Ion Conducting Polymer Electrolytes.

Single-ion conducting polymer electrolytes exhibit great potential for next-generation high-energy-density Li metal batteries, although the lack of sufficient molecular-scale insights into lithium transport mechanisms and reliable understanding of key correlations often limit the scope of modification and design of new materials. Moreover, the sensitivity to small variations of polymer chemical structures (e.g., selection of specific linkages or chemical groups) is often overlooked as potential design parameter. In this study, combined molecular dynamics simulations and experimental investigations reveal molecular-scale correlations among variations in polymer structures and Li+ transport capabilities. Based on polyamide-based single-ion conducting quasi-solid polymer electrolytes, it is demonstrated that small modifications of the polymer backbone significantly enhance the Li+ transport while governing the resulting membrane morphology. Based on the obtained insights, tailored materials with significantly improved ionic conductivity and excellent electrochemical performance are achieved and their applicability in LFP||Li and NMC||Li cells is successfully demonstrated.

Transport Properties of Waxy Crude Oil: A Molecular Dynamics Simulation Study.

To study the effects of paraffin on viscosity of waxy crude oil and transport properties of small molecules, light and waxy crude oil models were investigated at atmospheric pressure and 293-323 K temperature range using atomistic molecular dynamics simulations. The optimized parameters for liquid simulations all-atom (OPLS-AA) and atomistic polarizable potential for liquids, electrolytes, and polymers (APPLE&P) force fields were employed. The self-diffusion coefficients, viscosity, and paraffin configurations were compared for two oil models and between the two employed force fields. However, the behavior of paraffin molecules predicted by two force fields was quite different. Simulations using the OPLS-AA force field showed crystallization of longer paraffin molecules below 323 K, while simulations with the APPLE&P force field demonstrated a homogeneous mixture down to 293 K. To provide additional validation of the employed force fields, the density, diffusion coefficient, and crystallization of pure alkanes were compared with experimental data. The density and diffusion coefficients of n-C6 and n-C14 simulated with the APPLE&P force field were found to be in much closer agreement with the experimental data. The OPLS-AA force field was found to overestimate the crystallization temperature of pure alkanes. Therefore, simulations with the APPLE&P provide more realistic description of the waxy oil structure and transport properties. In this temperature range, the paraffin molecules are homogeneously distributed in the mixture, and viscosity of the system increased by a factor of two compared to light oil. Crystallization of paraffins requires lower temperatures or/and the presence of other components such as nanoparticles or asphaltene molecules to facilitate nucleation.

First-principles experimental demonstration of ferroelectricity in a thermotropic nematic liquid crystal: Polar domains and striking electro-optics.

We report the experimental determination of the structure and response to applied electric field of the lower-temperature nematic phase of the previously reported calamitic compound 4-[(4-nitrophenoxy)carbonyl]phenyl2,4-dimethoxybenzoate (RM734). We exploit its electro-optics to visualize the appearance, in the absence of applied field, of a permanent electric polarization density, manifested as a spontaneously broken symmetry in distinct domains of opposite polar orientation. Polarization reversal is mediated by field-induced domain wall movement, making this phase ferroelectric, a 3D uniaxial nematic having a spontaneous, reorientable polarization locally parallel to the director. This polarization density saturates at a low temperature value of ∼6 µC/cm2, the largest ever measured for a fluid or glassy material. This polarization is comparable to that of solid state ferroelectrics and is close to the average value obtained by assuming perfect, polar alignment of molecular dipoles in the nematic. We find a host of spectacular optical and hydrodynamic effects driven by ultralow applied field (E ∼ 1 V/cm), produced by the coupling of the large polarization to nematic birefringence and flow. Electrostatic self-interaction of the polarization charge renders the transition from the nematic phase mean field-like and weakly first order and controls the director field structure of the ferroelectric phase. Atomistic molecular dynamics simulation reveals short-range polar molecular interactions that favor ferroelectric ordering, including a tendency for head-to-tail association into polar, chain-like assemblies having polar lateral correlations. These results indicate a significant potential for transformative, new nematic physics, chemistry, and applications based on the enhanced understanding, development, and exploitation of molecular electrostatic interaction.

Molecular dynamics simulations of isothermal reactions in Al/Ni nanolaminates.

Molecular dynamics simulations of reactions in Al/Ni layered systems have been carried out under isothermal conditions for a wide range of temperatures and several system sizes. An embedded atom method potential, known to reasonably reproduce the phase behavior of Al/Ni, was employed. Simulations revealed reaction mechanisms involving an initial fast process and much slower more complex longer-time reactions. The initial reaction process consists of diffusion of Ni from the pure solid Ni phase into the molten Al phase, resulting in the formation of an Al-rich Al/Ni liquid. The initial reaction ends when the Al/Ni liquid becomes saturated in Ni and solid Al/Ni phases begin to form at the interfaces between the pure solid Ni phase and the Al/Ni liquid. The growth of these solid phases is intrinsically slow compared to the formation of the liquid and is further slowed by the need for Ni to diffuse through the growing interfacial Al/Ni solid phases. Analysis of the initial Al/Ni liquid forming process indicates Fickian behavior with the Ni diffusion coefficient exhibiting Arrhenius temperature dependence. The longer-time slow reaction process(es) resulting in the growth of Al/Ni solid phases do not lend themselves to detailed numerical analysis because of the complex dependence of the Ni transport on the detailed nature of the interfacial layers.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields.

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

How efficient is Li+ ion transport in solvate ionic liquids under anion-blocking conditions in a battery?

An experimental analysis based on very-low-frequency (VLF) impedance spectra and the Onsager reciprocal relations is combined with advanced analysis of dynamic correlations in atomistic molecular simulations in order to investigate Li+ transport in solvate ionic liquids (SILs). SILs comprised of an equimolar mixture of a lithium salt with glyme molecules are considered as a promising class of highly concentrated electrolytes for future Li-ion batteries. Both simulations and experiments on a prototypical Li-bis(trifluoromethanesulfonyl)imide (TFSI) salt/tetraglyme mixture show that while the ionic conductivity and the Li+ transport number are quite high, the Li+ transference number under 'anion-blocking conditions' is extremely low, making these electrolytes rather inefficient for battery applications. The contribution of cation-anion correlation to the total ionic conductivity has been extracted from both studies, revealing a highly positive contribution due to strongly anti-correlated cation-anion motions. Such cation-anion anti-correlations have also been found in standard ionic liquids and are a consequence of the constraint of momentum conservation. The molecular origin of low Li+ transference number and the influence of anti-correlated motions on Li+ transport efficiency have been investigated as a function of solvent composition. We demonstrate that Li+ transference number can be increased either by reducing the residence time between Li+ and solvent molecules or by adding excessive solvent molecules that are not complexing with Li+.

The 1-ethyl-3-methylimidazolium bis(trifluoro-methylsulfonyl)-imide ionic liquid nanodroplets on solid surfaces and in electric field: A molecular dynamics simulation study.

Atomistic molecular dynamics simulations were conducted to study the wetting states of 1-ethyl-3-methylimidazolium bis(trifluoro-methylsulfonyl)-imide ionic liquid (IL) nanodroplets on surfaces with different strengths of van der Waals (VDW) interactions and in the presence of an electric field. By adjusting the depth of Lennard-Jones potential, the van der Waals interaction between the solid surface and ionic liquid was systematically varied. The shape of the droplets was analyzed to extract the corresponding contact angle utilized to characterize wetting states of the nanodroplets. The explored range of surface-IL interactions allowed contact angles ranging from complete IL spreading on the surface to poor wettability. The effect of the external electrical field was explored by adding point charges to the surface atoms. Systems with two charge densities (±0.002 e/atom and ±0.004 e/atom) that correspond to 1.36 V/nm and 2.72 V/nm electric fields were investigated. Asymmetrical wetting states were observed for both cases. At 1.36 V/nm electric field, contributions of IL-surface VDW interactions and Coulombic interactions to the wetting state were competitive. At 2.72 V/nm field, electrostatic interactions dominate the interaction between the nanodroplet and surface, leading to enhanced wettability on all surfaces.

Charge Transport in [Li(tetraglyme)][bis(trifluoromethane) sulfonimide] Solvate Ionic Liquids: Insight from Molecular Dynamics Simulations.

Molecular dynamics simulations using fully atomistic polarizable force field have been performed on solvate ionic liquids (SILs), comprised of tetraglyme (G4) solvent molecules, Li+ cations, and bis(trifluoromethane) sulfonimide (TFSI) anions, [Li(G4)][TFSI]. The SILs with equimolar salt:G4 composition were investigated in the 303-373 K temperature range, whereas several systems with lower salt concentrations were investigated at 373 K. The simulations using polarizable force field demonstrate very good consistency of structural and dynamic properties with experimental data. The ability to accurately sample the ion transport mechanisms is particularly encouraging, taking into account that previous simulations employing nonpolarizable models had challenges in sampling dynamics in these systems. Here, we correlate Li+ ion local environment and glyme conformations with dynamic characteristics, such as residence time of species around Li+, self-diffusion coefficients, transference number, and conductivity. The analysis of contributions to Li+ mobility due to changing its local environment (i.e., moving from one glyme/anion to another) and from translational motion of Li+ with its' coordination environment showed significant dominance of the latter. The contributions of cross-ion dynamic correlations to the total conductivity have been quantified, showing strongly positive contribution from the cation-anion anticorrelation. Despite the high degree of Li-TFSI dissociation and positive contribution of the cation-anion anticorrelated motion to conductivity, the Li+ transference numbers for equimolar SILs are very low under the anion blocking conditions.

Role of cationic groups on structural and dynamical correlations in hydrated quaternary ammonium-functionalized poly(p-phenylene oxide)-based anion exchange membranes.

Extensive atomistic molecular dynamics (MD) simulations employing a polarizable force field have been conducted to study hydrated anion exchange membranes comprised of a poly(p-phenylene oxide) (PPO) homopolymer functionalized with quaternary ammonium cationic side groups and hydroxide anions. Representative membranes with different cationic structures have been investigated to study correlations between polymer architecture, morphology and transport properties of hydrated membranes. Specifically, hydrated polymers with five different quaternary ammonium cationic groups (R1: -CH3, R2: -C2H5, R3: -C3H7, R4: -C6H13 and R5: -C4H8OCH3) and degree of functionalization of 50% were investigated at three hydration levels (λ = Nwater/Ncation = 5, 10 and 17). Effects of the polymer structure on the distribution of water-rich domains and dynamic relaxations were systematically investigated to uncover the complex interplay between the degree of hydrophobicity/hydrophilicity of the cationic groups, morphology, connectivity of water domains, and the hydroxide transport mechanisms. Structural and dynamical analysis indicates that the bottlenecks, formed between the water-rich domains, create a substantial free energy barrier for hydroxide transport associated with the partial loss of anion hydration structure. The energy penalty associated with the loss of the hydration structure hinders the vehicular transport of the hydroxide anion. The optimal structure of functionalized homopolymer chains should be sufficiently hydrophobic to create nanophase segregation and form an interconnected network of water channels with a minimal amount of narrow bottlenecks that inhibit the vehicular motion of hydrated anions. We demonstrate that utilization of asymmetrically modified cationic groups is a promising route to achieve the desired water channel morphology at low hydration levels.

Grotthuss versus Vehicular Transport of Hydroxide in Anion-Exchange Membranes: Insight from Combined Reactive and Nonreactive Molecular Simulations.

Combined reactive and nonreactive polarizable molecular dynamics simulations were used to probe the transport mechanisms of hydroxide in hydrated anion-exchange membranes (AEMs) composed of poly(p-phenylene oxide) functionalized with the quaternary ammonium cationic groups. The direct mapping of membrane morphologies between two models allowed us to investigate the contributions of vehicular and Grotthuss mechanisms in hydroxide motion and correlate these mechanisms with the details of local structure. In AEMs with nonblocky polymer structure, where anion transport occurs through narrow (subnanometer size) percolating water channels, simulations indicate the importance of the Grotthuss mechanism. In nonreactive simulations, in order to diffuse through bottlenecks in the water channels, the hydroxide anion has to lose part of its hydration structure, therefore creating a large kinetic barrier for such events. However, when the Grotthuss mechanism is involved, the hydroxide transport through these bottlenecks can easily occur without loss of anion hydration structure and with a much lower barrier.

Application of Screening Functions as Cutoff-Based Alternatives to Ewald Summation in Molecular Dynamics Simulations Using Polarizable Force Fields.

The range-dependent screening of the charge-charge, charge-induced dipole, and induced dipole-induced dipole interactions was examined for a variety of liquids modeled using polarizable force fields. A cutoff-based method for calculation of the electrostatic interactions in molecular dynamics (MD) is presented as an alternative to Ewald-type summation for simulations of the disordered materials modeled using many-body polarizable force fields with permanent charges and induced point dipoles. The proposed approach was tested on bulk water, room-temperature ionic liquids, and solutions of ions in polar solvents. The smooth, short-range, and atom-type independent screening functions for interactions between the charges and induced dipoles were obtained using the force matching approach. An excellent agreement for both the magnitude and directionality of forces, structural and dynamic properties, was found in MD simulations utilizing the developed screening functions, compared to those with Ewald summation. While similar in shape and range, the charge-charge screening functions were somewhat dependent on the material chemistry. On the other hand, the charge-induced dipole and induced dipole-induced dipole screening functions were found to be nearly universal for the tested materials.

Multiscale Modeling of Structure, Transport and Reactivity in Alkaline Fuel Cell Membranes: Combined Coarse-Grained, Atomistic and Reactive Molecular Dynamics Simulations.

In this study, molecular dynamics (MD) simulations of hydrated anion-exchange membranes (AEMs), comprised of poly(p-phenylene oxide) (PPO) polymers functionalized with quaternary ammonium cationic groups, were conducted using multiscale coupling between three different models: a high-resolution coarse-grained (CG) model; Atomistic Polarizable Potential for Liquids, Electrolytes and Polymers (APPLE&P); and ReaxFF. The advantages and disadvantages of each model are summarized and compared. The proposed multiscale coupling utilizes the strength of each model and allows sampling of a broad spectrum of properties, which is not possible to sample using any of the single modeling techniques. Within the proposed combined approach, the equilibrium morphology of hydrated AEM was prepared using the CG model. Then, the morphology was mapped to the APPLE&P model from equilibrated CG configuration of the AEM. Simulations using atomistic non-reactive force field allowed sampling of local hydration structure of ionic groups, vehicular transport mechanism of anion and water, and structure equilibration of water channels in the membrane. Subsequently, atomistic AEM configuration was mapped to ReaxFF reactive model to investigate the Grotthuss mechanism in the hydroxide transport, as well as the AEM chemical stability and degradation mechanisms. The proposed multiscale and multiphysics modeling approach provides valuable input for the materials-by-design of novel polymeric structures for AEMs.

Supramolecular Self-Assembly of Methylated Rotaxanes for Solid Polymer Electrolyte Application.

Li+-conducting solid polymer electrolytes (SPEs) obtained from supramolecular self-assembly of trimethylated cyclodextrin (TMCD), poly(ethylene oxide) (PEO), and lithium salt are investigated for application in lithium-metal batteries (LMBs) and lithium-ion batteries (LIBs). The considered electrolytes comprise nanochannels for fast lithium-ion transport formed by CD threaded on PEO chains. It is demonstrated that tailored modification of CD beneficially influences the structure and transport properties of solid polymer electrolytes, thereby enabling their application in LMBs. Molecular dynamics (MD) simulation and experimental data reveal that modification of CDs shifts the steady state between lithium ions inside and outside the channels, in this way improving the achievable ionic conductivity. Notably, the designed SPEs facilitated galvanostatic cycling in LMBs at fast charging and discharging rates for more than 200 cycles and high Coulombic efficiency.

Structural and Dynamical Properties of Tetraalkylammonium Bromide Aqueous Solutions: A Molecular Dynamics Simulation Study Using a Polarizable Force Field.

Understanding the behavior of aqueous solutions containing tetraalkylammonium (TAA) cations is of great significance in a number of applications, including polymer membranes for fuel cells. In this work, a polarizable force field has been used to perform atomistic molecular dynamics (MD) simulations of aqueous solutions containing tetramethylammonium (TMA) or tetrabutylammonium (TBA) cations and Br counterions. Extensive MD simulations of TMA-Br/water and TBA-Br/water systems were conducted as a function of solution composition (ion pair:water molar ratios of 1:10, 1:20, 1:30, 1:63, and 1:500) at atmospheric pressure and 298 K. Our simulations demonstrate excellent agreement with available experimental data for solution densities and diffusion coefficients of different species as a function of solution composition, providing us confidence in analyzed structural and dynamic correlations. Various ion-ion and ion-water spatial distributions and the extent of cation aggregation are discussed in light of changes in the structure of cations hydration shells. The delicate balance between cation ionic core interactions with water and the hydrophobic interactions of alkyl tails leads to nontrivial self-assembly of TAA cations and the formation of an interpenetrating cationic network at higher concentrations. The ions and water dynamics are strongly coupled with the observed structural correlations and are analyzed in terms of various residence time, diffusion coefficients, and ionic conductivity.